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ABSTRACT 

A sample of ~1,800 halo subdwarf stars with radial velocities and proper motions is assembled from 
Bramich et al.'s (2007) light-motion catalog. This is based on the repeated multi-band Sloan Digital 
Sky Survey photometric measurements in Stripe 82. Our sample of halo subdwarfs is extracted via a 
reduced proper motion diagram and distances are obtained using photometric parallaxes, thus giving 
full phase space information. The tilt of the velocity ellipsoid with respect to the spherical polar 
coordinate system is computed and found to be consistent with zero for two of the three tilt angles, 
and very small for the third. We prove that if the inner halo is in a steady-state and the triaxial velocity 
ellipsoid is everywhere aligned in spherical polar coordinates, then the potential must be spherically 
symmetric. The detectable, but very mild, misalignment with spherical polars is consistent with the 
perturbative effects of the Galactic disk on a spherical dark halo. Banana orbits are generated at the 
1:1 resonance (in horizontal and vertical frequency) by the disk. They populate Galactic potentials 
at the typical radii of our subdwarf sample, along with the much more dominant short-axis tubes. 
However, on geometric grounds alone, the tilt cannot vanish for the banana orbits and this leads to 
a slight, but detectable, misalignment. We argue that the tilt of the stellar halo velocity ellipsoid 
therefore provides a hitherto largely neglected but important line of argument that the Milky Way's 
dark halo, which dominates the potential, must be nearly spherical. 

Subject headings: subdwarfs — Galaxy: kinematics and dynamics — Galaxy: structure — Galaxy: 

halo 



1. INTRODUCTION 

The kinematics of any stellar population are often most 
conveniently described by its velocity dispersion tensor 



(1) 



where the subscript indices denote one of the orthogo- 
nal coordinate directions, and the angled brackets rep- 
resent averag ing over the phase space d istribution func- 
tion (see e.g.. lBinnev fc Tremainel[2008D . The dispersion 
tensor is a symmetric second-rank tensor and so may al- 
ways be diagonalized. The principal axes of the tensor 
then form a velocity ellipsoid, which need not be aligned 
with the coordinate directi ons. However, as already re- 
ahzed bv lEddingtonl (|1915| ). if the gravity field is time- 
independent, then the alignment of the velocity ellipsoid 
is a powerful global probe of the gravitational potential. 
The triaxiality of the local halo v e locity ellipsoid 



is we ll -establ i shed (e .g.. IWoollevI Il978t IChiba fc Beers 
2005 iGouldl [2001 iKeolev et all l2007t iSmith et al " 



2003). However, despite its importance, the alignment 
of the velocity ellipsoid of halo stars has received very 
little attention. Here, we construct an unprcccdentcdly 
large sample of '^1,800 halo subdwarf stars with known 
distances, radial velocities and proper motions in Sloan 
Digital Sky Survey (SDSS) Stripe 82. We find that the 
velocity dispersion tensor is anisotropic, with an align- 



ment very close to that of the spherical polar coordinate 
axes. 

If the velocity ellipsoid is exactly aligned radially every- 
where, then the potential of the Milky Way is spherical 
in shape. The very mild misalignment that we detect is 
consistent with the influence of the Galactic disk on an 
underlying spherical Galactic halo potential. 

2. THE SUBDWARF CATALOG 

2.1. Sample Construction 

We construct our sample of subdwarfs us ing dat a from 
the s ixth SDSS data release (jAdelman-Mc Carth v et al.l 
l2008f). in particu l ar ut ilizing the light-motion catalog 
of iBramich et all (|2008[ ). This is built from the multi- 
epoch, multi-band {u, g,r,i, z) photometry available for 
one of the SDSS equatorial stripes (Stripe 82) and covers 
~ 250 deg^ in the right ascension range 2tf.^7 < a < 3^3 
and in the declination range |i5| < 1?26. A f ull description 
of ou r subdwarf sample is given elsewhere ([Smith et al. I 
[2009h . Here, we give a brief outline of the selection pro- 
cedure. 

The reduced proper motion is defined as 



5 logi 
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(2) 



Institute of 
bridge, Madingley 



Astronomy, University of Cam- 
Road, Cambridge CBS OHA, UK; 
msmithSast . cam . ac . uk , nweOast . cam . ac . uk 

^ Dark Cosmology Centre, Niels Bohr Institute, University of 
Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen 0, 
Denmark jinOdark-cosmology .dk 

^ Niels Bohr International Academy, Niels Bohr Institute, 
University of Copenhagen, Blegdamsvej 17, DK-2100 Copenhagen 
0, Denmark 



^mas yr 

where fi is the proper motion and r is the apparent 
magnitude i n the r band, correct ed for extinction using 
the maps of lSchlegel eFaTI ([19980 . The reduced proper 
motion is useful because it is independent of distance, 
namely 

Vtan 



Hr= Mr + 5 logio 



(3) 



.4.74 km s^i. 

where Mr is the absolute magnitude in the r band and 
^tan is the tangential component of the velocity with re- 
spect to the line-of-sight between the Sun and the star. 
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The reduced proper motion diagram is a plot of Hr ver- 
sus color g — i, in which the populations of disk dwarf 
stars, white d warfs and the halo subdwarfs are all sepa- 
rated (see e.g.. lVidrih et al.ll2007l . who have already con- 
structed a reduced proper motion diagram for Stripe 82 
to isolate ultracool and halo white dwarfs). 

Here, we are interested in selecting a clean sample of 
halo subdwarfs, and so we apply two pre-selection cuts 
in order to reduce the contamination.^ First, we only 
use stars that pass a quality cut such that uncertainties 
in the proper motion are smaller than 4 mas yr^^. Sec- 
ondly, the r magnitude of the star should be brighter 
than r = 19.5. This gives us a sample size of 372,811. 
The latter cut allows us to remove interlopers which may 
be at large distances and hence have small Hr despite 
having fj, ^ mas yr~^. Neither cut introduces any kine- 
matic bias. Note that we avoid cutting on ^ directly in 
order to simplify our calculation of the detection effi- 
ciency. 

In Figure [l] we show the reduced proper motion dia- 
gram. Owing to larger tangential velocities, the halo sub- 
dwarfs are clearly differentiated from the slower moving 
disk dwarfs. We estimate the location of the subdwarf 
boundary as 

Hr < 2.85(.g + 11.8 for {g - i) < 2 

Hr < 5.63(g + 6.24 for {g ~ i) > 2 

Hr > 2.85(.g + 15.0 for {g - i) < 1.3 

Hr > 5.63(g + 11.386 for [g - i) > 1.3, (4) 

and further clean the sample by rejecting all objects, 
which are within AHr = 0.5 of the boundary. These 
cuts result in a sample of 30,760 stars. We then cross- 
match with the out put from the SP SS SEGUE spectral 
parameter pipehne ()Lee et al.ll2008[ ). This provides us 
with the radial velocities for members of the subdwarf 
sample (with median error of 4.6 km s^^), although it 
significantly reduces the sample size to 2,210 stars. The 
SEGUE spectroscopic target selection is a complicated 
function. Suffice it to say that, despite its complexity, it 
is believed to be free of any significant kinematic biases. 

In order to recover the full six-dimensional phase-space 
informa tion, we mus t determine the distance to each sub- 
dwarf. Ilvezic etldl (|2008f l have already constructed a 
photometric parallax relation from SDSS observations 
of clusters, obtain i ng an intrinsic scatter of ^ . 2 ma g 
(jlvezic et al.l l2008t lJuric et all (20081 : ISesar et al.ll2008f ). 
Their parallax relation is also a function of [Fe/H], 
which we obtain from the spectral parameter pipeline 
(with median er ror of 0.14 dex). O verall, the errors ob- 
tained from the Ilvezic et al.l (|2008[ ) relation are impres- 
sively small, ai d ed by the high precision photometry from 
iBramich et al.l (|2008( ) - for our sample the median rela- 
tive distance error is 10.1%. 

In conjunct io n wit h the good accuracy on the 
IBramich et al] (|2008D proper motions (less than 
4 mas yr~^), we are able to construct an unprece- 
dentedly large sample of halo subdwarfs with accurate 
positions and kinematics. The median error on each of 
our velocity components is ~ 30-40 km s~^. Restricting 

Possible contaminants such as white dwarfs, d isk dwarfs or 
backg round giants can be considered to be negligible: ISmith et al. I 
((2003) conclude that the level of contamination for this sample is 
<1%- 
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Fig. 1. — Reduced proper motion diagram of Stripe 82, where 
the color-scale corresponds to number density (scaled so that the 
peak is unity). The solid lines show the location of our subdwarf 
boundary, while the dashed lines show the region adopted to re- 
duce contamination. To improve the clarity of the figure, we have 
incorporated a cut on the proper motion (p > 30 mas yr~^); the 
sample used throughout the paper has no such cut on /i. Note that 
the color-scale saturates at 20% of the peak density. 



ourselves to subdwarfs with heliocentric distances less 
than 5 kpc gives a final sample of 1,782. These stars 
lie at Galactocentric cylindrical polar radii between 7 
and 10 kpc, and at depths of 4.5 kpc or less below the 
Galactic plane. 

2.2. Kinematic Bias Quantification 

A drawback to our sample is that it is not kinemati- 
cally unbiased, as the cut in the reduced proper motion 
diagram is implicitly a function of the kinematics. There- 
fore, we have to model and understand the effect of this 
if we are to investigate the distributions of velocities in 
our sample. Notice that the kinematic bias comes solely 
from our cut on the reduced proper motion Hr - which 
actually selects stars via their tangential velocity rather 
than their proper motion. This makes the task of quan- 
tifying the bias significantly easier since we do not need 
to make any assumptions about the underlying distance 
distribution (i.e., luminosity function). 

We calculate our detection efficiency as follows. For 
each subdwarf in our final sample, we take the sky co- 
ordinates and create a mock sample of 50,000 fake stars. 
Then, for each mock star, we select Mr and (g — i) at ran- 
dom from our observed distributions. Note that for each 
realization, both the magnitude and color are assigned 
from one star, i.e., we do not assign an absolute magni- 
tude from one star and a color from another. We then se- 
lect kinematics for each mock star u sing the Ga lactoccn- 
tric halo velocity distrib utions fromlKeplev et a l. (200^ 
but with no net rotation (jAUende Prieto et al.l[2006l) . To 
calculate Vtan from the Galactocentric velocities requires 
us to assign a distance to each mock star, which we do at 



Tilt of the Halo Velocity Ellipsoid 



3 



random from the observed distribution. This means that 
the efficiency does have a dependence on the distance dis- 
tribution. However, this dependence is very mild since 
equation Q is a function of the tangential velocity rather 
than the proper motion. The efficiency is then given by 
the fraction of mock stars which pass our reduced proper 
motion cut. 

In order to check whether our results are dependent 
on the assumed halo velocity distri bution, we repeat th e 
calculations using the values from iKeplev et al.l ()2007D , 
but now assuming a rotational velocity of the halo of 
~ 20 km s~^ (in the direction of disk rotation). We find 
that uncertainties in the efficiency make little difference 
to the final determination of the tilts. 

3. THE ALIGNMENT OF THE HALO VELOCITY 
ELLIPSOID 

3.1. Method 

With the efficiencies in hand, we can now calculate the 
misalignment of the velocity ellipsoid of the SDSS subd- 
warfs. To do this, we transform our subdwarf velocities 
into Galactocentric spherical polar coordinates: is the 
radial velocity with respect to the center of the Galaxy, 
vg is the zenithal component measured from the North 
Galactic Pole, and Vfj, is the azimuthal component mea- 
sured such that the Galactic rotation has negative w^. 

The misalignment from the spherical polar coordinate 
surfaces can then be described by the correlation coeffi- 
cients and the tilt angles using the following formula 



where 



N 



and 



tan(2ai 



24 



(5) 



(6) 



Here the tilt angle corresponds to the angle between the 
i-axis and the major axis of the ellipse formed by pro- 
jecting the three dimensional velocity ellipsoid onto the 
ij-plane. (see e.g., Binney & Merrifield 1998, or Ap- 
pendix A of this paper) . We use the tilt angles to specify 
the orientation of the velocity ellipsoid, instead of alter- 
natives such as the Euler angles because the former are 
a natural extension of the familiar two dimensional case 
and much easier to visualize than other options. See also 
recent examples of using the tilt angles in similar context 
by Dehnen & Binney (1998) and Siebert et al. (2008). 

The measured sample covariance is due to both the 
true underlying covariance and the correlated measure- 
ment uncertainties, i.e., 



C0Vm[Vj,Uj] 



Cov[(5vi, 5vj 



(7) 



where Gov,, 



Vj\ is the covariance as measured from 
the sample and Cov[(5fi, bvj\ is the covariance of the error 
distributions. To account for the detection efficiency, we 
also calculate the sample covariances weighted by the 
inverse efficiency, e.g.. 



N 



G0Vni[w,;,Wj] = ■j^X!^'^ (^'^'^ ~ k^iM)) {Vj,k - {Vj,k)), 

(8) 



k=l 



fe=l 

and the summation is over the number of subdwarfs in 
our final sample. The weights Wk are proportional to 
reciprocal of the efficiency e^^ (normalized to unity 

such that X^feLi ''^k = !)■ 

As indicated by equation ([5]) , when the covariances are 
calculated we subtract the mean velocities, i.e., our re- 
sults are not dependent on the assumed refer e nce fr ame. 
However, as can be seen from iSmith et al.l (|2009t) our 
sample displays no significant net motion. 

3.2. The Tilt Angles and Correlations 

Let us first look at the tilt angles in the sample ignoring 
any stars with height \z\ < 1 kpc, so as to minimize 
(although not eliminate) any effect due to the Galactic 
disk potential. This gives us a sample of 1,532 subdwarfs. 
We find that the correlations and tilts are 



Coir[vr,vg] ^ 0.078±0.029, 
Gorr[z;^,W0]= -0.028±0.039, 
Cott[v^,V9]= -0.087 ±0M7, 



are = 3?4±1?3 

ar^ = -2?2±3?3 (9) 

ase = -37?4±20?4 



where the errors are obtained using the bootstrap tech- 
nique. We find no evidence of any clear tilt in the 
and a 09 terms. However, the tilt angle are is measured 
to be non-zero at about 3-(t level, though it is still very 
small. The good alignment of the velocity ellipsoid in 
spherical polars is also apparent from the velocity distri- 
butions in the {vr,ve) and {vr,v^) planes illustrated in 
Figure [21 in which the dashed lines show the orientation 
of the tilts. 

Note that the tilt angle a^e is not well constrained 
since the dispersions in the two components are similar 
and the covariance is small. Hence, the calculation of 
a^e involves the division of one small number by another 
small number, giving a large error. However, this term 
is not so important for the purpose of constraining the 
shape of the Galactic potential (see the next section). 

The area of sky covered by the Stripe 82 catalogue is 
249 deg^ . It is reasonable to ask whether the velocity dis- 
persion ellipsoid locally could have a different behaviour 
than that averaged over a larger region, especially as this 
effect has been observe d recently in numerical simula- 
tions (jZemp et al.ll2009t ). To study this, we investigate 
the variation in tilt as a function of height from the plane 
for \z\ < 4.5 kpc in Figure [3l Neither of the tilt angles 
are and show any significant trends. Although there 
are one or two isolated points at which the tilt angles 
lie at ~ 1-cr from the mean, they are localized in longi- 
tude, indicative probably of kinematic substructure and 
streams which are expected in the stellar halo in hierar- 
chical assembly models. The nature of the substructure 
in the SDSS subdw arfs is addressed in detail elsewhere 
(jSmith et al. Il2009[ ). We have also carried out this test 
by splitting our sample into three according to right as- 
cension, and find that the values of the tilt angles are 
consistent with those derived for the full sample. 

4. DISCUSSION 
4.1. Spherical Alignment 
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Fig. 2. — The efficiency corrected velocity distributions in the (vryVg) and {vr,v,p) planes for the sample of 1,532 subdwarfs with 
1 kpc < Izl < 4 kpc. The dashed lines show the orientation of the tilts. The apparent non-Gaussianity in the {vr,v,p) distribution is due 
to the variation of the efficiency correction across this plane; for certain regions of velocity space (particularly for fa —200 km s~^), our 
efficiency is low and hence the few stars that fall in this range are corrected by a large amount. This results in the narrow peaks in density 
in this figure, although when averaged over a laxger area the density is consistent with a Gaussian distribution. 




Fig. 3. — The variation of the tilt angles a^-g (black) and ct^^ 
(red) as a function of height z from the Galactic plane. The vertical 
error bars are l-tr, whilst the horizontal error bars give the bin 
width. Although there arc one or two locations where the deviation 
from the mean is greater than l-cr, there arc no obvious trends 
discernible. 



The ahgnment of the velocity elhpsoid of halo stars in 
spherical polar coordinate has substantial implications 
for the overall potential of the Galaxy. This constraint 
does not come from the Jeans equations, which merely re- 
quire that t he momentum flux ba l ances the gravitational 
forces (see lEvans. An fc Walked (|2009f) or IXn fc Evaiii 



(|2009( l for recent applications). Rather, the constraint 
comes from the deeper requirement that a phase space 
distribution function must exist. The theorem has been 
known for some time, although its wid espread app l icabil - 
ity has been obscure d by th e fact that lEddington ()1915D 
and IChandrasekhail (jl939[ ) introduced un necessary as- 
sumpt ions in its proof, as realized first by iLvnden-Belll 
(fT962 ). 

Let us first note that there are a number of trivial ways 
that permit the velocity dispersion tensor to be aligned 
in spherical polar coordinates. The simplest is to ask for 
the distribution function to depend on energy E alone, 
in which case the velocity dispersion tensor is everywhere 
isotropic. Or, we could insist that the distribution func- 
tion is given as / = f{E, \L\) in a spherical potential, or 
/ = f{E^ Lz) in an axisymmetric potential. Here, L is 
the angular momentum, whilst is the component of L 
that is parallel to the symmetry axis. All these options 
are not available to us because of the well-known and 
long established triaxiality of the velocity dispersion of 
halo stars with cr^^ > cr^^ > agg - see for example recent 
determ inations by iKeplev et al.l (|2007| ) or ISmith et al. I 
dlQOi). 

To generate the observed triaxial anisotropy of the 
velocity dispersion tensor, the phase space distribution 
function must depend on at least three independent in- 
tegrals of motion /.;, one of which may be chosen to be 
the same as the Hamiltonian, Ii = E{v'^,Vg,v'^;r,6,(f)). 
Here, it is assumed that the reference frame is chosen 
such that there is no net bulk streaming motion, that 
is, {vr) = {ve) = {v^) = 0. Therefore, cr^^ = (w^), 
a'^g ~ (vrVg) and so on. Accordingly, if the cross-terms 
(vrVg), and (vrV^) vanish everywhere, then any addi- 
tional integrals of motion that isolate the distribution 
function must be even functions of Vr- This is equiva- 
lent to the statement that the u^-dependence of the in- 



Tilt of the Halo Velocity Ellipsoid 



5 



tegrals is only through the square of the radial veloc- 
ity component w^, and independent of the sign of Vr- 
Then, the second isolating integral I2 can be recast as a 
globally defined function independent of using the en- 
ergy integral, that is I2 ~ l2{E;vg,Vcj,;r,0,(j)). However, 
since I2 is an integral of motion, its Poisson bracket with 
the Hamiltonian must vanish, from which it follows that 
I2 al so must be independent of r as well (jLvnden-Belll 
Il962[ ) . Hence, the integral can always be cast in the form 
I2 = l2iE;vg,v^;6,4i) and so is completely independent 
of both r and (its conjugate momentum) Vr- 

This implies that the radial coordinate in the 
Hamilton-Jacobi equation must separate and so the po- 
tential has to be of the form 

^(r,0,0)=V^oM + ^M, (10) 

where ipQ and '5 are arbitrary functions of the indicated 
arguments. In fact, by exactly the same line of reasoning 
applied to the third isolating integral I3, if (vgv^) also 
vanishes, then the potential has to have the form 

^(r, e, </>) =Mr) + ^ + (11) 



7 sm U 

although this is a stronger result than we will need here. 

If the potential has the form (|10p . then Poisson's equa- 
tion implies that the total density of stars and dark mat- 
ter is 

p(r,0,0)=po(r) + ?M. (12) 



where 



1 



Po 



AttGt'^ dr 



dr 



1 d 



sm( 



(9* 



1 52* 



That is to say, the dipole potential is necessarily asso- 
ciated with an astrophysically unrealistic density cusp 
diverging as unless $7 = 0. Consequently, we have 
= 0.^ This leaves us with the theorem that: 

// a steady state stellar population has a 
non-degenerate (i.e. triaxial) velocity dis- 
persion tensor whose eigenvectors are every- 
where aligned in spherical polar coordinates, 
then the underlying gravitational potential 
must he spherically symmetric. 

This theorem is implicit in iLvnden-Belll (IT96^, wh ereas 
restricted versions we re known to lEddingtonl ()1915f ) and 
IChandrasekhail (| 19391 ). In fact, from the preceding proof, 
the crux of the result lies in the existence of the second 
integral that forces the separation of the radial part of 
the Hamilton-Jacobi equation, and the presence of the 
third integral is of a secondary importance. Hence, the 
theorem can be relaxed somewhat: 

// the potential is non-singular, it is a suffi- 
cient condition for a spherical symmetry that 
one of the non- degenerate eigenvectors of the 
velocity dispersion tensor is aligned radially 
everywhere. 

^ Mathematically, 't may be any combination of spherical har- 
monics Y^{d,(f>) with / = 1. However, there must not be any 
physical source for the pure-dipole gravitational field, and so ^ 
may be set to zero by appropriate choice of the coordinate origin. 



This follows, as spherical symmetry is guaranteed for 
non-singular potentials once the radial part of the 
Hamilton-Jacobi equation can be separated off. 

Of course, spherical alignment of the velocity disper- 
sion tensor does not imply that the stellar density it- 
self is spherically symmetric. In fact, many investigators 
have already found th at a^g ^ a'^^ for samples of halo 
stars fe.g.. lGould|[2003l : iKeplev et al1l2007l : ISmith et al.1 
l2009( l. which implies that the stellar halo has a den- 
sity distribution that is flattened or triaxial even though 
the gravity field is nearly spherical. T his is the case , 
for example, in th e models presented bylWhitd ([l985h . 
lArnoldl (|1990D . and lEvans. Hafner fc de Zeeuwl (|1997t ). in 
which the phase space distribution function depends on 
energy and the all three components of angular momcn- 
tum, / = /(£;,£). 

4.2. Nearly Spherical Alignment 

The tilt angles and a^g arc consistent with zero, 
but our results give a very small, but non-zero measure- 
ment of the tilt angle a^e (at the S-cr level) . In an exactly 
spherical potential, the (vrVg) cross-term and hence are 
must vanish. However, even if the dark halo is spher- 
ical, the Galactic potential is not spherical due to the 
influence of the bulge and the disk. Of course, at the 
distances probed by our sample of SDSS subdwarfs, it 
is the halo that dominates the gravitational potential, 
whilst the disk gives the main perturbation. 

The main consequence of the disk is to convert 
the planar rosette orbits of a spherical potential into 
the short-axis tube orbits of a mildly oblate poten- 
tial. However, in typical axisymmetric Galactic po- 
tentials, the key 1:1 resonance (in cylindrical coordi- 
nates R and z) also occurs at Galactocentric radii be- 
tween 2 and 10 kpc, causing the axial orbits to be- 
come unstable to out-of-planc perturbations, and sir- 
ing the family of banana or saucer orbits (see e.g. , 
Pfennig er 1984; Miralda-Escude fc Schwarzschildl Il989l: 
Schwarzschildl Il993t lEvand Il994l : iBinnev fc Tremaind 
20081 Section 3.7.3). It is to these two orbital fami- 
lies - the short-axis tubes and the bananas - that our 
SDSS subdwarfs will belong. 

We can confir m this by computing o rbits in the Milky 
Way potential of lFellhauer et all (j2006[ ) , which comprises 
a spherical isothermal halo and a Miyamoto-Nagai disk. 
Our initial conditions were chosen as follows: since the 
potential is axisymmetric, we start our orbits off from a 
location a; = 8 kpc, y ~ kpc and with z distributed 
uniformly from —1 to —4 kpc (which is a reasonable ap- 
proximation of our observed distribution). Our initial 
velocities were chosen according to the trivariate ellip- 
soidal Gaussian ar = 142 km s^^, = 81 km s^^, and 
ag ~ 77 km s^^, i.e., according to the valu es derived for 
our ob served sample of halo subdwarfs by ISmith et al. I 
(|2009f) . This is of course not a bona fide distribution 
function, but just a convenient sampling function to scan 
phase space. The orbits were calculated using a fourth 
order Runge-Kutta method with a time-step of 0.1 Myr, 
which results in a fractional energy change of better than 
10^^ over a typical orbital period of ^ 200 Myr. The 
orbits so produced are indeed predominantly short-axis 
tubes (~ 80%), but with a reasonable mixture of banana 
orbits (~ 20%). Some examples arc illustrated in the 
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Fig. 4. — Cross-sections in the meridional plane of typical orbits of the SDSS subdwarfs; leftmost two panels showing thin and thick tube 
orbits, rightmost two panels showing banana or saucer orbits. 



panels of Figure |4l 

It is immediately apparent on geometrical grounds that 
the banana orbits can only give positive contributions to 
(vrVg). By contrast, the short-axis tubes can yield both 
positive and negative contributions. Hence, the effect of 
the Galactic disk is to create a family of orbits for which 
the tilt term a^g cannot exactly vanish. 

In fact, it is easy to bolster this geometric argument 
with a quantitative calculation. Using our sampling func- 
tion, we generated 96 initial conditions and integrated 
the orbits for 1 Tyr, recording the velocities each time 
a test particle returned to its initial position (defined 
in practice to within a spherical volume of 300 pc ra- 
dius). The enormous timescale is to maximize the num- 
ber of times a particle returns to its initial position, so 
that a^g for the orbit can be calculated as accurately as 
possible. Then, we average the resulting values of a^g 
for every crossing and for all of our orbits. This gives 
(vrVg) = 690 km^ s~^, which is very comparable to our 
observed covariance of 893 ± 335 km^ s~^. 

However, if we divide our orbits into banana and short- 
axis tubes (banana orbits are defined as those which are 
asymmetric about z = or the Galactic plane) and recal- 
culate the average tr^g, then we find that the 19 banana 
orbits have (vrVg) = 1740 km^ s~^, while the remaining 
77 short-axis tube orbits have {vrVg) = —140 km^ s~^. 
The banana orbits indeed contribute almost all of the 
effect. As a check, if we repeat this procedure in a spher- 
ical potential (i.e., without the Galactic disk), we indeed 
find that (vrVg) is consistent with zero, as it should be. 

A very mild misalignment in the velocity ellipsoid of 
the SDSS subdwarfs with respect to spherical polar co- 
ordinates is therefore naturally explained by the effects 
of the Galactic disk. 

5. CONCLUSIONS 

The Milky Way's dark halo dominates the gravity field, 
and different datasets have yielded very different values 
for its flattening. For example, analyses of the varia- 
tion of the thickness of the Galaxy's gas layer with ra- 
dius point to an oblate halo with axis ratio q ~ 0.7 
()OUing fc Merrifieldl[200l[ l. This is consistent with the 
typical values of halo flattening measured for some ex- 
ternal galaxies, using the dynamics of polar rings or the 
shapes of the isophotes of the X- ray emission for early- 
type galaxies. On the other hand. lFtellhauer et al.l ()2006l ) 
reached the very different conclusion that the dark halo 
of our own Galaxy must be very close to spherical, based 



on dynamical modelin g of th e bifur cation in the Sagit- 
tarius stream. Finallv. fHelmil ()2004l ) used the kinematic 
data of stars in the leading arm of the Sagittarius to 
argue that the dark halo is prolate. 

Here, we have drawn attention to another - somewhat 
neglected - probe of the shape of the dark halo, namely 
the tilt of the velocity ellipsoid of halo stars. In a pure 
spherical halo, the velocity ellipsoid of any tracer popu- 
lation ~ regardless of whether its density distribution is 
flattened or triaxial - must be aligned in the spherical po- 
lar coordinate system. By contrast, the behavior of the 
velocity ellip soid in flattened halos is muc h more varied. 
For example. iLevison fc Richstond (|l985al lbl) constructed 
a number of halo models with ellipticities between 0.3 
and 0.7 using the Schwarzschild method. They show nu- 
merous plots of the behavior of the velocity ellipsoid, 
from which it is apparent that it deviates, often very 
strongly, from spherical polar alignment. 

In this paper, the tilt of the velocity ellipsoid of the 
Milky Way halo stars has been measured using an un- 
precedentedly large sample of '^1,800 halo subdwarfs ex- 
tracted from Bramich et al.'s (2008) light-motion cata- 
log, itself based on the repeated Sloan Digital Sky Survey 
observations of Stripe 82. We find that the velocity el- 
lipsoid of the halo subdwarfs is very closely aligned with 
the spherical polar coordinate system. In particular, two 
of the tilt angles are consistent with zero, whereas the 
final tilt angle ~ corresponding to the (vrVg) term - is 
very small. We have shown that this effect is consistent 
with the influence of the Galactic disk, which must cause 
a mild asphericity in the total potential even if the dark 
halo itself is spherical. This asphericity gives rise to a 
family of orbits, the banana or saucer orbits, which are 
typically represented in samples of stars at the locations 
of the SDSS subdwarfs. On geometric grounds alone, 
(vrVg) does not exactly vanish for banana orbits, and so 
their small admixture causes a very mild misalignment. 

Our argument here is open to the criticism that we 
have not shown that the velocity dispersion tensor is 
aligned everywhere in spherical polar coordinates. We 
have merely shown that the alignment is very close to 
spherical polars for halo subdwarfs at heliocentric dis- 
tances of < 5 kpc along the ~ 250 deg^ covered by SDSS 
Stripe 82. Nonetheless, this is still a striking and un- 
expected result over a swathe of Galactic locations that 
provides a new line of attack on the awkward question 
of the shape of the Milky Way's dark halo. It would 
be very interesting to extend our results to samples of 
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halo stars in different directions. This is particular the 
case for high latitude samples of halo stars, for which 
any effects due to the Galactic disk are negligible, and so 
for which alignment with the spherical polar coordinate 
system should be perfect. 
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APPENDIX 

THE TILT ANGLES AND THE VELOCITY ELLIPSOID 
After the mean motion has been subtracted, the velocity dispersion is 

4= Jd^vv.v.f. (Al) 

If we consider a coordinate transformation u^' = Ai'jVj (henceforth the summation convention for repeated indices is 
used) in velocity space given by a matrix in SO (3), it is straightforward to verify that the velocity dispersion transforms 
as a component of a tensor, i.e.. 



2 



where the Jacobian relating the coordinate transformation is unity. 

Note from the definition that the velocity dispersion tensor is symmetric and its three eigenvalues are all positive 
definite. Hence, it is also invertible. Let {i:fj) = {<jfj)~^ be the matrix inverse of the velocity dispersion tensor. The 
velocity ellipsoid is then defined to be a quadric surface in velocity space 

<;fjViVj = 1. (A3) 

Next, from the rule involving the inversion of the product of matrices and the fact A^/j is an orthogonal matrix (i.e., 
A~j} = Aju), it is clear that (;fj also transforms as a tensor, that is to say, (,f,j, = Ai/feAj/;<j^;. Then, 
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and thus the equation of the velocity ellipsoid is form-invariant under orthogonal coordinate transformations. If 
the primed coordinate axes are chosen to be aligned to the principal axes of the velocity ellipsoid, the equation of 
the ellipsoid in the primed coordinate is in the canonical form. It then follows that ^.f,^, and subsequently af,ji are 
diagonalized in the same primed coordinate system. In other words, the coordinate transformation to the system 
with the basis set given by the principal axes of the velocity ellipsoid is identical to the one diagonalizing the velocity 
dispersion tensor. 

Next, we consider the projection of the velocity ellipsoid. Here, we suppose that the projection is onto the v^Vy- 
plane, although the other two follow basically the same procedure. Let us think of a line with fixed {v^, Vy) but with 
varying v^- If {vx,Vy) falls within the projection of the ellipsoid, the line intersects the ellipsoid in two points. On 
the other hand, if the point (vx,Vy) is outside the projection in the u^Uy-plane, the line does not cross the ellipsoid. 
Following similar logic, we find that the point {vx, Vy) is on the projection of the ellipsoid in the VxVy-pleLnc if the line 
running perpendicular to the plane through the given point is a tangent to the ellipsoid. Now, if the equation of the 
ellipsoid is considered as a quadratic equation for z given fixed {vx,Vy), the preceding argument indicates that setting 
its discriminant to be zero defines the equation of the projection in the VxVy-plane. Therefore, we find the equation of 
the projection of the velocity ellipsoid onto the WajWy-plane; 

Cyyvl-2CxyVxVy + Cxxvl=';zz^ (A4) 

where Cxx, Cxy, and Cyy are the matrix cofactors of (r^^,, <j^y, and c;^^, respectively. This traces an ellipse in VxVy-p\au.e. 

The projected velocity dispersion tensor corresponding to this ellipse may be defined analogously to the 3-d case. 
After some algebra, we find that 

(^yyl'^zz '^^xy/'^zz \ _ -'^ f ^xx Cxy 1 _ f '^xx ^xy 1 ( A K\ 

-Cxyl^lz Gxxl^lz J ~ d\ Cxy J " [ 0% oly J ^^^^ 

where D — |^,?^ | is the matrix determinant of iSij)- In other words, the projection of the velocity ellipsoid onto the 
Uj^Uy-plane is the same as the 2-d velocity 'ellipsoid' calculated with the 2x2 {vx^ 7;j,)-submatrix of the original 3x3 
velocity dispersion matrix. Consequently, the tilt angle defined as in equation ([6]) is the same as the angle between 
the principal axis of the ellipse that is the projection of the velocity ellipsoid onto the ij-plane and the coordinate 
axis. We note however that this is not the same as the angle between the coordinate axis and the projection of the 
principal axis of the velocity ellipsoid onto the same plane - that is to say, the principal axis of the velocity ellipsoid 
is not necessarily projected into the principal axis of the projected ellipse. 



